Entanglement, number fluctuations and optimized interferometric phase measurement 
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We derive a phase-entanglement criterion for two bosonic modes which is immune to number fluc- 
tuations, using the generalized Moore-Penrose inverse to normalize the phase-quadrature operator. 
We also obtain a phase-squeezing criterion that is immune to number fluctuations using similar 
techniques. These are utilized to obtain an operational definition of relative phase-measurement 
sensitivity, via analysis of phase measurement in interferometry. We show that these measures are 
proportional to enhanced phase-measurement sensitivity. The phase-entanglement criterion is a 
hallmark for a new type of quantum squeezing, namely planar quantum squeezing. This has the 
property that it squeezes two orthogonal spin directions simultaneously, which is possible owing 
to the fact that the SU(2) group that describes spin symmetry has a three-dimensional parameter 
space, of higher dimension than the group for photonic quadratures. The practical advantage of 
planar quantum squeezing is that, unlike conventional spin-squeezing, it allows noise reduction over 
all phase-angles simultaneously. The application of this type of squeezing is to quantum measure- 
ment of an unknown phase. We show that a completely unknown phase requires two orthogonal 
measurements, and that with planar quantum squeezing it is possible to reduce the measurement 
uncertainty independently of the unknown phase value. This is a different type of squeezing to the 
usual spin-squeezing interferometric criterion, which is only applicable when the measured phase is 
already known to a good approximation, or can be measured iteratively. As an example, we calcu- 
late the phase-entanglement of the ground state of a two-well, coupled Bose-Einstein condensate, 
similar to recent experiments. This system demonstrates planar squeezing in both the attractive 
and repulsive interaction regimes. 



I. INTRODUCTION 

Entanglement criteria are widely used to identify non- 
classical resources for potential applications in quantum 
technology. One application is the enhancement of mea- 
surement sensitivity. In practice, the most sensitive mea- 
surements are often interferometric. Hence, the measure- 
ment of an unknown quantity is reduced to the mea- 
surement of a phase-shift. In this communication, we 
analyze how non-classical, entangled states can increase 
phase-measurement sensitivity. To achieve this, we will 
introduce both an operational measure of relative phase, 
and a corresponding signature of phase entanglement be- 
tween two Bose fields, using well-defined interferometric 
particle-counting procedures. This is shown to quanti- 
tatively measure the enhancement of an interferometric 
measurement. It is the interferometric equivalent of the 
spin-squeezing criterion [IH3 , which is known to measure 
the non-classical precision of a clock Q . 

We introduce a relative phase operator which is well- 
defined in the case of variable total particle number, by 
using the generalized inverse method to prevent singu- 
larities in the inverse number operator. It is the sim- 
plest relative phase operator that is measurable inter- 
ferometrically. Number fluctuations are always found 
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in current experimental photonic and atom interferom- 
eter phase measurements. Hence, we clarify the opera- 
tional phase-measurement procedure already used heuris- 
tically to analyze experiments 0-0|- Most importantly, 
we show that when the measured phase is unknown 
prior to measurement, a state preparation that involves 
mode-entanglement is optimal, and is closely related to 
the relative-phase operator. This complements previous 
studies which generally assume either that the phase is 
already known to a good approximation, or that the 
phase-shift remains constant during repeated measure- 
ments. Here, we use the minimal number of measure- 
ments possible. This is inherently different to strategies 
employed to estimate phase through sequential or mul- 
tiple measurements, which assume the phase-shift is a 
classical, time-invariant quantity. 

Our entanglement measure is a normalized form of the 
recently introduced Hillery-Zubairy (HZ) non-hermitian 
operator product criterion Q , similar to that introduced 
in a previous paper We prove that this normal- 

ized form is a phase-entanglement signature for two Bose 
fields, and has the advantage of being almost immune 
to total number fluctuations. We show how this criterion 
can be interpreted as a variance measure which quantifies 
entanglement, and has a direct physical interpretation as 
the enhancement of phase measurement sensitivity in an 
interferometer. This is directly related to the idea of pla- 
nar quantum squeezing, in which the quantum noise is 
simultaneously reduced in two orthogonal directions of 
phase- measurement . 
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The present analysis focuses on linear, two-mode inter- 
ferometry with particle-counting detectors. This tech- 
nique is by far the most commonly used technique for 
phase- measurement, and therefore deserving of a care- 
ful analysis. Our results, while less general than an 
analysis of completely arbitrary quantum measurements, 
are able to be implemented immediately, since two-mode 
interferometers are widely available both for photonic 
and atomic fields. We discuss techniques for generating 
the required entangled input fields using atomic sources. 
As an example, we consider how to obtain these types 
of quantum state through the creation of a correlated 
ground-state in a coupled, two-mode Bose-Einstein con- 
densate (BEC) with either attractive or repulsive S-wave 
scattering interactions. We show that as well as giving 
sub-shot noise (squeezed) phase noise in two orthogonal 
phase directions simultaneously, it is possible to obtain 
nearly Heisenberg- limited peformance in one of two phase 
directions, which is useful when the phase is known ap- 
proximately. 



II. QUANTUM PHASE MEASUREMENTS 

Interferometers are designed to measure a relative 
phase-shift 0, typically between two beams. The phase 
information must first be encoded on quantum fields be- 
fore it is measured. Hence, there is a close relationship 
between interferometry, which measures a phase-shift in 
a medium, and the measurement of phase of a quantum 
field, which is where the interferometric phase informa- 
tion is stored. The relationship is more important than 
meets the eye. An experimentally measured phase-shift 
is can neither be truly classical nor time-invariant, which 
are commonly used assumptions. For a variety of rea- 
sons, repeated measurement is not always possible, and 
one must regard interferometric measurement as primar- 
ily a quantum measurement problem. 

In view of this, we first review earlier approaches to 
phase measurement of quantum fields, in order to explain 
and motivate the approach used in the rest of the paper. 



A. Phase operators 

We start with a generic problem in quantum phase 
measurement: how can one measure the phase of a quan- 
tum field or mode (a). Classically, one divides up a field 
amplitude into intensity (N a ) and phase (</> a ) by intro- 
ducing: 

a=y/N^e<*: (1) 

Next, if one measures the complex amplitude a, one sim- 
ply classically normalizes to obtain the classical phase, 
4> a = — i In \a/y/Na\ . We note that this gedanken- 
experiment for classical phase measurement involves 
measuring both real and imaginary components of the 
amplitude. 



Quantum studies of this generic problem date back to 
the early attempt of Dirac [lOj to define a quantum phase 
operator from a canonical commutation relation of form 



<l>a,N a 



-i. The underlying mathematical problem is 



that a hermitian quantum operator <f> a for the phase of 
a single harmonic oscillator operator a strictly does not 
exist, which has been the topic of many previous studies 
[lll - [l3j . This can be seem most easily from the commu- 
tation relations, [aa^ — a^d] = 1. If phase is hermitian, 



then d = V NU, where U is a unitary operator such that 
Wlf = UU^ = 1, and N = a^a. From the commutators, 
N - 1 = UNU^. This violates unitarity, since a unitary 
transformation cannot change an operator's eigenvalues. 

This problem does not of course occur classically, and 
it is at the root of quantum phase-measurement prob- 
lems. However, the idea of a phase-operator with canon- 
ical commutators is approximately valid at large particle 
number, and suggests the existence of a fundamental un- 
certainty principle called the Heisenberg limit: 

Acj)AN > -. 

Given that AN < N max /2, this has led to the idea of a 
Heisenberg limit of A<p > 1/N max on phase-measurement 
with at most N max particles. Sometimes one inter- 
changes TV and N max , as they are often related. 



B. Truncated Hilbert space methods 

One resolution of the lack of a hermitian phase opera- 
tor due to Pegg and Barnett, is to truncate the Hilbert 
space to a maximum boson number of s [T3 |. Next, one 
defines a phase eigenstate \6) as a discrete Fourier trans- 
form of number states \n), using: 



1 s 



e m6 In) 



(2) 



From this starting point, it is clear that a hermitian 
phase operator can simply be obtained from the defini- 
tion: 



a. = E 



'm | u ml \ w m \ 



(3) 



Here 8 m = #o + 27rm/ (s + 1), and 8o is a reference phase. 
This is a mathematically consistent approach, which re- 
solves the issue of hermiticity, but it leaves a number of 
practical questions unanswered. 

In particular, what is the physical meaning of the max- 
imum number s and reference phase 9^1 Is it possible to 
take the limit of large s in a unique way? From an op- 
erational perspective, what (if any) is the relationship 
between the abstract operator <j) a and an interferomet- 
ric measurement? One of the purposes of this paper is 
to understand how these questions can be answered and 
implemented using interferometry. 
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C. Relative phase operator 

Another solution along these directions is to define a 
relative phase operator, <fi = (p a ~<Pb, for two modes a and 
b [HI. This has the conceptual advantage that it corre- 
sponds to operational procedures - which always involve 
relative phase measurement - and clarifies the meaning of 
the reference phase. A relative phase operator is consis- 
tent philosophically with the fundamental idea of relative 
measurements in physics, and is the most operationally 
meaningful way to define quantum phase in many cases. 

With this approach, one works in a space whose al- 
gebra is defined by the equivalent angular momentum 
operators in the Schwinger representation: 



j 



i(atS-aSt), 



i 



J z = -{ a^a - tfb 



N = tfa + tfb. 



(4) 



It is usual to assume that one has an eigenstate of 
the total number N, with eigenvalue N and hence an 
equivalent angular momentum eigenstate of J — N/2. 
As we show later, this assumption is questionable in real 
experimental measurements. The states of well-defined 
phase then correspond to linear combinations of angular 
momentum eigenstates with J Z \J, m) — m\J, m), and 
one then has a physical phase basis of: 



\0) p 



1 



V2J 



1^ 



J, m) 



(5) 



m=—J 



This allows the definition of a relative phase- difference 
operator which is Hermitian, and has a discrete spec- 
trum. The fact that the angular momentum Hilbert 
space is finite provides a natural explanation of the trun- 
cation parameter s in the Pegg-Barnett approach. How- 
ever, it is still not immediately clear how to relate this ab- 
stract proposal to interferometric phase-measurements. 

We also note that in practical interferometry experi- 
ments, it is nearly impossible have an input state that 
has a well-defined total particle number, especially at 
large mean particle number. Instead, the most common 
situation is that there is an initial mixture of particle 
numbers, with number fluctuations that are typically at 
least Poissonian. We will return to the issues of interfer- 
ometric measurement and number fluctuations in later 
sections. 



D. Quantum sine and cosine operators 

An alternative resolution of the phase-measurement 
question is to define quantum sine and cosine operators. 
This is a way to reach the phase through measurement 



of the real or imaginary part of the amplitude, an idea 
which has a clear analog in the classical world, as de- 
scribed above. Operationally, this is the quantum ver- 
sion of a proposal by Zernike fl6j to analyze coherence 
in classical interferometry. The original idea of Zernike 
was to relate coherence properties directly to measured 
classical fringe visibility. 

Extendin g th is to the quantum theory of a single quan- 
tized mode [ll|, one can define a normalized amplitude: 



E 



1 + N 



-1/2 



(6) 



from which the sine and cosine operators are obtained 
via: 



*4 



E-& 



(7) 



This is still operationally unclear for general applications, 
except in the limit where one arm of an interferometer 
is a large 'classical-like' local oscillator. In this limit, 
the approach is closely related to the theory of optical 
quadratures, where the two operators involved have com- 
mutators similar to the quantum position and momentum 
operators: for a field mode o, the quadrature amplitudes 
are defined X = (a + a t )/2, and P = (a- a t )/ (2i). 

A drawback of the local oscillator approach implicit 
in this method is that in practical terms it is highly 
resource-hungry. A classical local oscillator must nec- 
essarily involve a very large boson number in the local 
oscillator mode, even though this mode is not formally 
included as part of the measured system. We note that 
there is a generic issue, which is that since these opera- 
tors do not commute, it is not possible to measure both 
quadratures simultaneously. 



E. Relative phase quadrature operators 

Alternatively, the same general idea of quadrature 
measurement can be applied to the relative phase be- 
tween two modes. This also allows a better understand- 
ing of the actual resources involved in the measurement, 
since it effectively combines both the measured beam and 
the local oscillator into the measured operator. The com- 
bined approach is studied in an early review of Carruthers 
and Nieto [l3[ , and reduces to defining the cosine and sine 
operators of the phase-difference as: 



Cl2 — C1C2 + Sx02i 
C\2 = S1C2 — S2C1. 



(8) 



This has the virtue that these operators commute with 
the total particle number, N = Ni + N2. This number is 
kept finite, and provides an indication of the true resource 
needed for the measurement. A related definition for a 
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Figure 1. Interferometric measurements can be reduced in the 
simplest case to a measurement of outputs from a beamsplitter. 



two-mode BEC was proposed by Leggett [l7|, who sug- 
gested defining a phase operator following the approach 
of Carruthers and Nieto, except that: 



N/2 - J ; 



(N/2- 



.V 



-1/2 



X 



(9) 

Here the phase is encoded in the rotations in the X — Y 
plane, a convention we will follow in this paper unless 
otherwise noted. In these two approaches, the proposed 
quadrature operators were not analyzed from the per- 
spective of interferometric fringes and their noise proper- 
ties. We will show that the operationally most relevant 
approach to interferometry measurements is to use a dif- 
ferent normalization to either of the above suggestions. 



III. PHASE OPERATORS AND 
INTERFEROMETRY 

An operational analysis of interferometric measure- 
ments can be reduced in the simplest case to a measure- 
ment of outputs from the final beam-splitter as shown in 
figure [T] with one mode b experiencing an the unknown 
phase-shift of 4>, while the other mode a is shifted by a 
fixed reference phase 9. The other components of the in- 
terferometer are then part of the quantum state prepara- 
tion that determines the output expectation values. This 
is shown schematically as the input beam-splitter with in- 
put modes ai, bi, in the schematic Mach-Zehnder (MZ) 
diagram in figure [TJ 

In the following, we shall mostly focus on the gen- 
eral scheme in which the intermediate modes a, b have 
an arbitrary state preparation. However, we shall also 
treat particular examples where the state preparation is 
obtained through the Mach-Zehnder protocol. In this 
case we analyze the state preparation of modes ai, bi, 
as a practical route towards preparing the intermediate 
modes, and also introduce an additional phase-shift on 
the input to the MZ, for reasons explained in the last 
section. 



A. Quantum limits to classical phase-estimation 

The problem of measuring a quantum phase <f> is re- 
lated, but not identical, to the analysis of quantum lim- 
its to estimation of a classical phase. There is an es- 
sential difference, since classical phase-estimation usually 



assumes there is a phase-shifting element that produces 
a well-defined classical phase-shift c\>, which is supposed 
to be time- independent. Many treatments assume that 
the phase can be measured iteratively without disturb- 
ing it, to improve the accuracy. This limitation rules 
out many situations where the phase evolves in time - or 
where the phase experiences a back-action which changes 
the phase after measurement. Other approaches to the 
problem make the assumption that it is possible to con- 
struct arbitrary quantum states and measuring devices. 
As a result, these treatments generally are not applica- 
ble to two-mode interferometry, although they may be 
applicable to some future phase-measuring device. 

We first introduce earlier approaches to phase- 
estimation. Pioneering work by Caves [Hj], who treated 
bosons in the context of gravity-wave detection, showed 
that two-mode interferometry sensitivity could be im- 
proved above the shot-noise level i.e., the Standard Quan- 
tum Limit (SQL) of A</5 = 1/VN: 



A</5 < !/■< 



(10) 



This required non-classical, 'squeezed state' input radia- 
tion, reaching a maximum sensitivity near the Heisenberg 
limit of 



A<f> = l/N 



(11) 



for an input state with N average particle number. 
Squeezed states allow the uncertainty of one observ- 
able to be reduced below the standard quantum limit, 
at the expense of the complementary observables, so 
that the Heisenberg uncertainty relation is still satisfied 
[l9j . Thus, for single mode optical amplitudes where 
X = (a + at)/2, and P = (a - &)/(2i), for which 
AXAP > 1/4, squeezing of X occurs when AX < 1/2. 
This is clearly very closely related to the quadrature 
phase-operator [Tl| approach of equation (0. 

Paradoxically, the usual squeezed state technologies of 
parametric down-conversion are rather inefficient. The 
total resources employed, in terms of boson number prior 
to down-conversion, are generally no better than coher- 
ent interferometry [20j for a given phase sensitivity. As 
pointed out by Caves, there is still an advantage of this 
method for gravity-wave detectors, as it reduces the back- 
action caused by radiation pressure. 

The treatment of Caves also assumed prior knowledge 
of the unknown phase, to allow a linearized treatment in 
the limit of small fluctuations. This type of analysis was 
applied to fermion interferometry by Yurke I2ll. |22l . and 
was later extended to multiple measurements [23j, with 
similar conclusions, except for the replacement of N by 
the total number of particles involved, Ntot- 

Squeezed quantum fluctuations have been shown to 
enhance the sensitivity of other sorts of measurements 
[H, [2f|. Wineland et. al., in the context of atomic 
clock measurements Q, showed that there was a close 
relationship between interferometry and the concept of 
spin-squeezing [l|. The uncertainty relation for spin is 
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AJ Y AJ Z > |(j x )|/2, and spin squeezing exists when 
the variance of one spin is reduced below the standard 
quantum limit (SQL) 



AJ Y < J\(J x )\/2. 



The spin squeezing factor 



2JAJ Y 
\(J X )\ 



(12) 



(13) 



was introduced for a collection of N two-level atoms or 
equivalently, for two occupied modes for which a col- 
lective pseudo-spin is defined, and J = N/2. In spec- 
troscopy or interferometry, the final measurement is that 
of a spin component, measured as a population difference 
between the two quantum states. 

Since we wish to consider phase in the initial state 
prior to a beamsplitter, we consider that the x direction is 
chosen to be that of the large spin vector, so that ( J x ) ~ 
J. Fluctuations are thus spin-squeezed when A 2 J Y < 
J/2 = N/A. The precision of the quantum measurement 
is given by 0, [26[ 



A(t> = ts/VN (14) 
which is enhanced over the SQL (fTUl) when 

& < 1 (15) 

and reaches the Heisenberg limit as £s — > 1/ \Z~N. This 
implies that there are large fluctuations in J z . The most 
extreme case of this is when we choose the input state 
to be the eigenstate of the two-mode phase operator, Eq 
([5|), which establishes a connection between these two 
approaches. In this case, one finds that 



2 jZ 



A* J 



N 
12 



[N + 2] 



(16) 



Further analysis has suggested ways to reach these 
Heisenberg limits (fTTj) through macroscopic superposi- 
tions or 'NOON' states (26[. More recently, nonlinear in- 
terferometry [27j was suggested as a route to go beyond 
the Heisenberg limit, although this requires specific non- 
linear couplings. We note that all of these techniques use 
a linearized approach to phase-estimation. This means 
that the phase must already be known to an excellent 
approximation prior to measurement. In practice, this 
implies repeated measurements to refine the estimation 
of the phase until the final, high-precision measurement 
is made. In particular, we show below that a two-mode 
phase eigenstate only gives a low interferometric mea- 
surement variance when the phase is known almost per- 
fectly prior to the measurement. 

There are general treatments that typically involve 
more than single-pass, two- mode interferometry [28[. 
Sanders and Milburn [29j found the optimal measurement 
and state to determine the phase <f>, based on the two- 
mode phase-operator approach of equation ([5]). Berry 



and Wiseman [30i . |31| showed that this canonical mea- 
surement cannot be realized by counting particles in an 
interferometer (figure [IJ , and proposed alternative iter- 
ative schemes. Such idealized measurements have also 
been analyzed by many authors using techniques like the 
quantum Fisher information (32l . [33| . 

The assumption of a time-invariant, classical phase- 
shift used in iterative schemes can rule out such tech- 
niques for many applications. Phases often vary in time, 
and one cannot always avoid quantum back-action. The 
physical reason for this is that a phase-shift corresponds 
to an energy-shift term in the radiation-matter Hamilto- 
nian [34| . This therefore has dynamical consequences for 
the object being measured, and may change the phase. 
Hence, the use of iterative or repeated measurement 
schemes is not always feasible. Strategies of this type 
are different to the present theory, which focuses on min- 
imal numbers of measurements using a two-mode linear 
interferometer. 



B. Phase- measurement operator 

We now return to the fundamental question of how 
to measure a phase-shift using a quantum phase- 
measurement operator. We wish to use a strictly oper- 
ational definition, solely utilizing the interferometer out- 
puts. If we define the operator phases so that c,d = 
\ae~ 10 ± be~ l 0~\ /y/2, then the measured outputs of the 
quantum interferometer in terms of a, & are: 



N+ 



c f c ± dU 

1 -N± 1 - 

2 2 



ae 



(17) 



Comparing these quantities with the equivalent angular 
momentum operator approach from equation (j4]), we see 
that these quantities can be rewritten as: N± = N/2 ± 
Jy, where N is the total number operator, and: 

J = J x cos(0 -9) + J Y sin(0 - 6). 

In any given measurement, there are always two outputs 
which can be measured simultaneously. The quantum 
phase of the output field - which is also related to the 
unknown phase-shift </> - can be estimated from the ratio: 



N N 



J (<f>) = lim ■ 
£^0+ 



2 ( e + TV 2 



(18) 



Here N + is the Moore-Penrose generalized inverse [35| 
of N. This is a well-defined hermitian observable that 
commutes with J^, and gives the least-squares solution 
of any inversion or variational problem involving JV. The 
generalized inverse N + has many of the properties of a 
standard inverse, including the property that its eigenval- 
ues are iV _1 for number states with total number N > 0. 
However, it is zero (not infinity) for the vacuum state. 
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This means that, unlike the standard inverse, it has a 
well-defined value for all quantum states. More details 
are in the Appendix. 

In experiments designed for accurate phase measure- 
ment for large average particle numbers, events with zero 
total particle number occur with vanishingly small prob- 
ability, and generally can be neglected, in which case 
N+ w N^ 1 . Of course, this is only true approximately, as 
the number operator has no standard inverse. The quan- 
tity J (</>) is always measurable and hermitian. Hence, 
J (0) can be called a relative phase quadrature operator, 
and is the fundamental quantity measured in any phase- 
sensitive interferometric experiment. 

It is vital to normalize by the particle number at each 
measurement - as indicated in the above operator - for 
the simple reason that in general, the particle number 
is not known in advance. It is often theoretically as- 
sumed that the total particle number is known prior to 
measurement. This is rarely found in real experiments, 
especially as the number is increased. An inspection of 
the experimental protocols actually used in recent BEC 
interferometry experiments [H, [H, 0] shows that the op- 
erator given above corresponds rather closely to the way 
that data is analyzed in practice. Our analysis therefore 
provides a theoretical justification for these operational 
procedures. We will treat the effects of typical Poissonian 
particle number fluctuations in a later section. 

The operator J ((f)) is different to both the complex 
phase- difference quadrature operators of Nieto and Car- 
ruthers [lH, and of Leggett, owing to the type of nor- 
malization chosen here. If we compare the current ap- 
proach of equation (TT5)) with these earlier suggestions in 
equations ((8|) and ((9]), there is a very important differ- 
ence. The operator J (0) given above is uniquely defined 
for all inputs, and can be measured completely from a 
single, combined measurement of the two interferometer 
outputs. It is not obvious how one can measure the ear- 
lier proposed phase-measurement operators in practical 
interferometry experiments, since they appear to require 
the simultaneous measurement of non-commuting output 
operators like J x and J z . Of course, this does not rule 
out more sophisticated operational measurements, as the 
combined operators are hermitian; but these measure- 
ments do not appear feasible with simple beam-splitters 
and photodetectors. 

If multiple measurements are made sequentially, then 
more sophisticated iterative phase-estimation techniques 
are possible [H, [28|, [3(| [H, |37j]. As pointed out above, 
this does not help in experiments where the phase is 
changing. Often, only a single measurement is possi- 
ble. We also recall that many previous analyses are 
conditioned on having a priori approximate knowledge 
of the unknown phase-shift (f>. In the following, we fo- 
cus instead on optimizing the sensitivity of the opera- 
tional phase measurement equation (|18[) . for a range of 
unknown phase-shifts. In other words, we assume that 
the phase is known to lie in a given interval that is not 
vanishingly small. 



C. Entanglement and squeezing 

Interferometric sensitivity and particle entanglement 
have previously been linked through criteria involving 
Fisher information [3^]. Sorenson et. al. [38] have 
shown that a measure of particle entanglement is the 
spin squeezing criterion (|1 5|) : a fixed number N of two- 
level systems (spin- 1/2 particles) are separable when 
p = J2r-PrPi---Pn ar± d hence entangled when 



o < - 



NAJ Z / Y 



\(J X )\ 



< 1. 



(19) 



Here we define the collective spins associated with N spin 

1/2 systems: J 8 = J2k=i (° = X > Y > Z )> and J® is the 
spin of the fc-th particle. The Heisenberg Uncertainty 
Principle places a lower bound on £s, because of the fi- 
nite size of the Hilbert space. The precise values for 
the lower bound for fixed N have been determined by 
Sorenson and Molmer in Ref. (3^, and decrease with in- 
creasing spin J. The spin squeezing criterion has been 
measured experimentally in BEC interferometry [343; 
and is related to phase-measurement efficiency when the 
phase value is approximately known in advance Q, as 
summarized by (|14[) . 

In this paper, we will generalize this criterion to include 
number fluctuations, and also treat a very different type 
of entanglement, that between two distinct spatially sep- 
arated locations, rather than between many qubits. We 
will relate the sensitivity of the phase measurement (fT5| 
of an unknown phase to a special type of two-mode en- 
tanglement between a and b. Two-mode entanglement is 
defined as a failure of the separable model 



P : 



PRPaPb 



R 



(20) 



where Pr > 0, ^2r Pr = 1 and p^ b are density operators 
for states at a/b. 

Many criteria for two- mode entanglement exist, but we 
are interested only in those interferometric measures of 
entanglement that can enhance the phase measurement 
task of Figure Q] |H E3, E3 ■ The observab les that can 
be measured are given in equation (QJ and (|18l) . These 
expressions are written in terms of the modes a and b that 
are the inputs to the final beam splitter. One may also 
consider a MZ type of experiment with a phase-rotation 
and two beam-splitters as is illustrated in Figure [1] In 
this case, the MZ internal modes are related to the input 
modes a,, 6, by a = (— ia, + 6i)/v / 2, b = —(iai + 6i)/v / 2, 
and the measured outputs are simply rotated versions of 
the Mach-Zelmder input operators, with: 



J x 



J, 



z 

x x . 



(21) 



We see from these expressions ^ and (fT5)) that the 
sensitivity of the phase measurement will depend on the 



7 



noise levels of the two orthogonal components, j x N + 
and J N + (in terms of the interferometer modes a and 
b). It makes sense to then choose an input state for the 
interferometer that will maximally reduce the noise in 
both of these components simultaneously. In fact this 
requirement is closely related to an entanglement mea- 
sure. Hillery and Zubairy (HZ) showed [8j that for any 
separable state (|20p . 



A 2jx +A 2jY > ^ N y 2 . 



(22) 



Entanglement between modes a and b is thus detected 
when (El) fails: 



< E HZ = 



A 2 J X +A 2 ,P 

(N)/2 



< 1. 



(23) 



Ehz = 1 gives the Standard Quantum Limit (SQL) noise 
level, which is the noise level Ehz obtained when the 
modes a and b are in the separable product of coherent 
states, |a)|/3). It is not possible however to choose a state 
so that both variances A 2 J X , A 2 J Y are zero. 



D. Planar Quantum Squeezing 

If we consider the Heisenberg uncertainty principle in 
the X — Y plane, we see that it has the form A J Y A J x > 
\{J z )\/2. Here the optimal situation is obtained for 
equal beam intensities entering the beam-splitter, so that 
( J z ) = 0. This appears to provide no lower bound to the 
measured quadrature variances, and hence to the phase 
noise. However, appearances can be very misleading. In 
fact, a non-zero lower bound to Ehz exists, because the 
variances of J x and J Y cannot be simultaneously zero. 

For fixed N — 2 J, this bound has been determined. It 
is known that 



Cj/J < E HZ , 



(24) 



where the coefficients Cj ~ 3 (2 J) 2 ^ 3 /8 as J — ► oo [43 _ 
This means, however, that both the orthogonal variances 
in a phase-measurement can be simultaneously reduced 
below the shot-noise level, since we are minimizing the 
sum of the phase variances, A 2 J X + A 2 J Y . In gen- 
eral, a noise reduction of the sum of two variances below 
the shot-noise level is called planar quantum squeezing 
or PQS [43| . as it minimizes quantum noise in a plane, 
rather than just in one direction on the Bloch sphere. It 
has the advantage that noise reduction for phase mea- 
surement occurs regardless of the value of the unknown 
phase. 

It is instructive to compare the optimal PQS state with 
the relative phase eigenstate. The minimal variance PQS 
state is: 



i 



j 

V2TT1 ^ 7 

m— — J 



I J, m) 



(25) 



The asymptotic limit of the optimal coefficient R m , 
which minimizes the sum of the quadrature variances, is 
then a Gaussian of form: 



,-m 2 /{1u m ) 



(26) 



where the variance in the space of j z eigenvalues is 
a m = A 2 J Z = (j 2 /2) 2 ^ 3 . The other important prop- 
erties are that, to leading order: (J x ) ~ J,A 2 J X ~ 
(2 J) 2/3 /8, A 2 J Y ~ (2 J) 2/3 /4 . Thus, the optimal PQS 
state reduces the variance in J x and J Y simultaneously, 
with both variances well below the shot-noise level. 

By contrast, for a relative phase state at large J, 
we have ( J x ) ~ ttJ/4, A 2 J x ~ (2/3 - tt 2 /16)J 2 , 
A 2 J Y ~ yir7(J)> A 2 J Z ~ J 2 /3. Hence, the Heisen- 
berg uncertainty principle in the Z-Y plane is obeyed, 
since AJ Y AJ Z > \(J x }\/2. However, quantum squeez- 
ing only occurs in the Y spin direction, while in both 
the other spin directions the noise is greatly increased 
above the shot-noise level. This means that in an inter- 
ferometric measurement using a relative phase state, the 
reference phase-offset 9 must be adjusted to match <p with 
high precision, even though <f> is of course unknown prior 
to measurement. This adjustment is necessary to avoid 
contamination of the results with high levels of noise from 
the X spin direction, which are well above the shot-noise 
or Poissonian level. The underlying cause is that inter- 
ferometric measurements do not simply project out the 
relative phase eigenstates. 

We see that the main advantage of PQS states in in- 
terferometry is that it is possible to have sub-shot pre- 
cision in both the measured spin directions simultane- 
ously. This is advantageous when the measured phase 
is truly unknown. At first, it may seem that this is less 
than optimal as a squeezing strategy when the phase is 
known approximately. For the optimal PQS state de- 
scribed above, which minimizes the variance sum, neither 
of the variances are close to the Heisenberg limit. Impor- 
tantly, there are a range of possible PQS states in which 
the relative variances in the X and Y directions can be 
adapted to the desired measurement strategy, including 
states in which PQS - with the variance sum below the 
shot-noise level - is combined with nearly Heisenberg- 
limited variance reduction in one of the two directions. 
This possibility is discussed in the last section, together 
with practical techniques for achieving it. 



IV. OPERATIONAL CRITERIA AND NUMBER 
FLUCTUATIONS 

In practical interferometry, the total number of input 
bosons usually changes at each measurement. Hence the 
ensemble used for averaging has a finite distribution over 
the particle number. This is caused by a number of fac- 
tors. In optical lasers, it is caused by technical noise 
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in the optical pumping process, as well as well-known 
quantum noise effects during stimulated emission and 
out-coupling [Hj]. In BEC and atom lasers, the factors 
involved range from fluctuations in the initial atomic den- 
sity distribution in the magneto-optical trap, to quantum 
noise due to the atomic collisions that occur in the evap- 
orative cooling process [45j . These number fluctuations 
are due to the non-equilibrium mechanisms that gener- 
ate a laser or BEC respectively, and there is no reason to 
assume either a canonical or a grand canonical ensemble. 

The direct use of the HZ spatial entanglement crite- 
rion is highly sensitive to total number fluctuations. For 
this reason, we will introduce entanglement and spin- 
squeezing definitions that are normalized by the total 
number. We refer to such normalized entanglement mea- 
sures as phase- entanglement and phase-squeezing mea- 
sures, as they measure correlated and reduced noise 
phases. There is another possible strategy, which is to 
simply reject all measurements that have the 'wrong' par- 
ticle number. This allows a conditional number state 
measurement to be obtained. While this is feasible, it is 
also extremely inefficient, since most attempted measure- 
ments yield no information at all about the phase with 
this strategy. 



A. Experimental number nuctuations 

In the case where the number fluctuations are Pois- 
sonian, the probability that there are exactly N bosons 



to the effect of total number fluctuations. This is mis- 
leading: total number fluctuations do not destroy entan- 
glement. Accordingly, it is important to use entangle- 
ment and phase-sensitivity measures that allow for num- 
ber nuctuations. Our general criteria therefore include 
number fluctuations with an arbitrary variance. These 
criteria can still be used to describe idealized experiments 
without number fluctuations, even though this is not very 
realistic. In the next section, we will make use of a Pois- 
sonian distribution to model the behavior of typical BEC 
experiments, with a low atom number of N ~ 100. 



B. Entanglement and squeezing criterion 

To treat phase measurement including total number 
fluctuations, we have introduced normalized spin opera- 
tors: J = J 6 N + . Here, N + is the Moore- Penrose gen- 
eralized inverse of the number operator. Detailed prop- 
erties and proofs are given in the Appendix. We can now 
use these normalized operators to derive general opera- 
tional criteria for entanglement and squeezing, which ex- 
tend the results obtained above to a realistic environment 
with number-fluctuations. We note that somewhat dif- 
ferent results have been obtained previously with number 
fluctuations, in work that used un-normalized operators 

El. 



1. Phase-entanglement criterion 



1 



P{N) = — (N) N e 



NJ-(N)) 



(27) 



While this may not be the best input state for a phase 
measurement, this distribution does give a number stan- 
dard deviation of = \/iV", which is a typical or- 
der of magnitude for the number fluctuations in a well- 
stabilized photonic laser or BEC. This fluctuation is in 
the total number TV, prior to any beam-splitter. Interfer- 
ometric beam-splitters introduce relative number fluctu- 
ations in addition to the total number fluctuation. 

Highly-stabilized semiconductor lasers have reached 
slightly lower number variances than Poissonian, in a re- 
stricted frequency range |46j]. For an atomic BEC exper- 
iment, atom number statistics are difficult to measure 
accurately to this level of precision for large N. Num- 
ber fluctuations of at least the Poissonian level are found 
in almost all current BEC experiments [f| where data 
is available. Just as with lasers, it is possible to obtain 
lower variances than this, with some restrictions. In the 
best results to date, standard deviations as low as 0.6y/~N 
(below the Poissonian level) were observed at very small 
atom numbers of N « 60. Super-Poissonian variances 
were found for larger numbers of N > 500 (47| . 

When calculating the entanglement parameter Ehz 
for typical states with Poissonian fluctuations, we find 
that the apparent entanglement is greatly decreased due 



We now introduce a phase-entanglement measure that 
is robust against number nuctuations. We show in the 
Appendix that in a number-fluctuating environment, en- 
tanglement between modes a and b is confirmed via a 
phase-entanglement criterion that uses the generalized 
Moore- Penrose inverse of the number operator, N + : 



2 t y 



E 



i>ii 



A 2 J A + A 2 J 



< 1. 



(28) 



2. Phase-squeezing criterion 

Similarly, we show in the Appendix that entanglement 
between N spin- 1/2 systems is confirmed by a normalized 
spin squeezing criterion, which we term phase-squeezing: 



f Z/Y 



jZ.Y 



\(J X )\ 



< 1. 



(29) 



The two criteria (|2"5)) and (f2"9"|) and their application to 
determine the enhanced sensitivity of a two-mode atom 
interferometer, in particular a BEC atom interferome- 
ter for which incoming number fluctuations are included, 
form the major results of this paper. 
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C. Phase sensitivity 

Next, we will obtain a detailed understanding of the 
relationship between our phase-entanglement measure, 
and phase-measurement sensitivity. The crucial issue in 
phase measurement is the measurement sensitivity, or 
smallest measurable phase-shift. This is related to the 
differential signal to noise ratio, given by pH l22l | 



d 4 - m- 1 



f 



,d(J). 



(Ajy 



(30) 



The smallest measurable change in phase in a single mea- 
surement is A(f>. Figure[T]depicts an unknown phase shift 
(fi (to be measured) relative to a fixed phase shift 9. We 
suppose for simplicity that (a^6iV + ) = \(a^bN + )\, so that 
the direction of the mean spin of the state to be used in 
the interferometer will be along the X axis: i.e. when 
cb = 0: (J X N+) = \{tfbN+)\, (J Y N+) = (J Z N+) = 0. 
For a controlled reference phase shift 9' , two successive 
orthogonal measurement settings, 9' and 9' + tt/2, will 
allow a determination of the unknown phase (fi: 



(%>,0')) =cos(0 
(J(0,0'+7r/2))=-sin( 



- 9')\(a r bN+)\, 
- 9')\(a r bN+)\. 



(31) 



A single measurement setting 9' cannot determine the 
unknown phase completely, since the information given 
is about cos tp only. The mean differential signal for mea- 
surement j(<fi, 9') = J r >'N + is -(J X N+) siny, and A</> as 
given by ([50)1 for this measurement is (<p = <fi — 6') 

{A(fif = {(AJ X ) 2 cot 2 M + (AJ^) 2 }/|(J^)| 2 ,(32) 

together with a similar expression obtained in the orthog- 
onal direction. The objective is to determine the condi- 
tions on the interferometric state so that the uncertainty 
in the phase estimation is minimized. 

The Standard Quantum Limit (SQL) sensitivity A(f> = 
1 / vN, as given by equation (fT0[) , is obtained when fields 
a and b are formed from a number state |iV) incident at 
one port of a beam splitter . with a vacuum state input 
at the second port (Bl. l2ll. l4l|. An entangled state results 

[HI, for which (J x ) = N/2, (AJ X ) 2 = and {AJ Y f = 
N/4. Then, for all phases ip, it is readily shown that 
Acfi ~ 1/VN. 

For some entangled states, it is well-known jiol [50| | 
that the phase sensitivity can be enhanced below the 
SQL. The most well-studied cases however consider a 
small phase shift about a fixed phase reference [3J, |41| . 
It is evident from (|32|) that the maximum differential for 
j((f>,0') is at ip = tt/2, for which (A0) w/2 = AJ Y /\(J X )\. 
The sensitivity at this point is thus given by the nor- 
malized spin squeezing parameter (|29[) which reduces to 
(|19p for fixed number inputs. Sub-shot noise sensitiv- 
ity is achieved when A<fi <l/y (N), so by the definition 



(I2"9")l , sub-shot noise enhancement occurs for interferomet- 
ric states satisfying 



< 1. 



(33) 



The technique relies on an accurate estimate (fix of the 
unknown phase, combined with setting 9 to (fix — tt/2, 
so that subsequent measurements measure small shifts 
near the optimal phase. We will show in the next section 
that a near Heisenberg-limited sensitivity of (A0) ff /2 ~ 
0{s/2/N) is predicted for this case, when the two-modes 
a and b of an atom interferometer are prepared from a 
two-mode double-potential well BEC ground state. 



D. Estimation of an unknown phase 

The question of phase estimation with an unknown 
phase and limited number of measurements is a differ- 
ent issue [1^1 . Where we restrict to phase estimation via 
the interferometric scheme figure Q] based on the number 
difference measurements (|17M18[) , we see from (I3"2l that a 
noise-reduction enhancement over a range of angles, with 
a reduced variance in both AJ X and A J Y is needed. This 
is essential where there is no prior knowledge of the phase 
(fi and successive adaptive phase measurements [30j are 
not possible. We note that the sensitivity of the measure- 
ment J(<fi,9') is destroyed by the divergent contribution 
evident in (|3"2"j) at <p ~ 0, tt, unless A J x = 0, which 
places a severe limit on the interferometric state. This is 
not a necessary consideration, however, if the full phase 
is to be measured via both orthogonal measurements set- 
tings given by (f3"Tj) . The settings in the "quiet" quadrants 
6' = (((> + tt/2) ± tt/4 and 9' = (<f> - tt/2) ± tt/4 have 
enhanced sensitivity over those in the "noisy" quadrants 
9' = (fi±Tr/4 and 9' — (</>+7r)±7r/4, and for any unknown 
phase cfi one of the orthogonal settings 9' or 9' + tt/2 will 
be in the useful quadrants. Least squares estimation is 
an obvious strategy here. 

We thus consider the following strategy. The first read- 
ing of the pair J{<fi, 0) or J [cfi, tt/2) determines values for 
cos (fi and sin (fi, and thus the location of the phase in the 
plane. In this way, it can be determined which of J(<fi, 0) 
or J{(fi, tt/2) has measured in the quiet quadrants. Sub- 
shot noise sensitivity is then guaranteed at all unknown 
angles (fi, for this preferred measurement, provided it can 

be shown that Aifi < 1/ y (N) across the entire two quiet 
quadrants. According to (|32p . the worse-case sensitivity 
for these quiet quadrants is at <p — ±7r/4, ±37r/4, and is 



2 

max 



given by (A(fi) 
condition for (A<fi) 



A 2 (J X ) + A 2 (J Y ) /\(J X )\ . The 
to be sub-shot noise is (A(fi) w < 

1 / \J (N) which is quantified by a phase-sensitivity mea- 
sure: 



2 

max 



J (N)(N+)E ph /2 

Vph= — m\ — <L 



(34) 
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When the interferometric fields satisfy (f3~4")h sub-shot 
noise sensitivity for all angles <f> is guaranteed, for any 
fixed measurement setting 9 within the two quiet quad- 
rants for measurement of </>. 

The fundamental quantum limit for (|34j) is given by 
the smallness of r) p h, which is linked to the uncertainty 
relation for the sum of the two spins J x and J Y . It is 
therefore important to determine a tight lower bound on 
this sum, in order to obtain the ultimate phase interfer- 
ometric sensitivity. The real question becomes to what 
extent can we still minimize A 2 J Y , given that the sum 
A 2 (J X ) + A 2 (J Y ) is also to be minimized. The answer is 
not the same as for two complementary observables like 
two optical quadratures, or position and momentum, for 
which the commutator is a constant. The uncertainty 
relation for spin operators has a state-dependent form, 

AJ X AJ Y > \ (J z )\/2. 

which means the two variances AJ X and AJ Y can both 
be reduced below the shot noise level. For (J z ) = 
the Heisenberg uncertainty principle is unable to give 
any bound at all. This is possibly misleading, since the 
Heisenberg uncertainty principle is simply a bound that 
does not guarantee it can be saturated. As discussed in 
section MIDI a recent analysis of the spin- variance un- 
certainty leads to a tight bound on the variances: 

A 2 J X + A 2 J Y > Cj, 

where Cj is a function of the total spin J with an 
asymptotic limit of 3(2J) 2 / 3 /8. For J = N/2 one finds 
Acj> > y^UJ/J -> 0{VT!>/N 2 / 3 ), which we will show in 
the next section is predicted for the two- well BEC ground 
state, in the both attractive and repulsive regimes. The 
fundamental limit is below the SQL of 0(l/y/~N) over all 
the phase angles in the half-plane, but can only reach the 
Heisenberg limit of 0{\/N) over part of the range. We 
next present the details of how these levels of sensitivity 
can be realized in a BEC interferometer. 



V. BEC INTERFEROMETER 

We consider how the criteria (J25J), and ||31P for en- 
tanglement and phase-measurement sensitivity are satis- 
fied in a typical cold-atom experiment. The mechanism 
for two-mode squeezing and entanglement here is similar 
to that first realized for optical modes using four-wave 
mixing [5l| , except employing the ground-state of an in- 
teracting BEC. We consider an idealized two-mode or 
two-well BEC with a normalized self-interaction coeffi- 
cient g, and a linear tunneling coupling n between two 
modes with boson operators a and b, described by the 
Hamiltonian [52J 

H = n{a)b + a6 f ) + | [aVaa + tftfbb]. (35) 



Solutions for the two-mode entanglement of the ground 
state have been presented in Refs. [H, H^, for the case 
of an initial state of N atoms distributed evenly between 
the wells. 

In order to interpret these figures it should be kept 
in mind that the attractive BEC case generates a nearly 
perfect planar quantum-squeezed (PQS) state. This is 
close to an optimal minimization of the variance sum, 
reducing phase-noise in both quadratures. The resulting 
state is used in direct interferometry, with a single beam- 
splitter. On the other hand, while the repulsive BEC case 
also produces a PQS state, it has a different characteristic 
with one of the quadratures having a phase-noise level 
reduced almost to the Heisenberg limit, while the other 
phase quadrature is not at the Heisenberg limit. This 
state is used in Mach-Zehnder interferometry, with an 
additional phase rotation to generate the state with the 
optimal characteristics for phase- measurement. 



A. Entanglement criteria 

We present solutions for the ground state of ([331) . in- 
cluding number fluctuations as described in section HV Al 
The dashed curves in figure [5] show that the normalized 
phase-entanglement criteria Eq. (|25|) E p h < 1 detects the 
two-mode entanglement of the ground state of a two- well 
BEC [4]-|6|, l9( in a way that is almost immune to Poisso- 
nian number fluctuations. On the other hand, the solid 
curves for Ehz reveal that the entanglement detected via 
the un-normalized Hillery-Zubairy criterion Ehz < 1 is 
very easily destroyed by number fluctuations. 

The figure [3] shows that the normalized spin squeezing 
parameter also detects squeezing and particle entangle- 
ment in a way that is insensitive to number fluctuations. 
The solid curves plot the squeezing predicted for a fixed 
N = 100, where we recall that this parameter is not de- 
fined except in the case of fixed N, but that according 
to Eq. (fTi?|) will detect both squeezing and entanglement 
among the N particles when £g < 1. The dashed curves 
plot the results for squeezing of the normalized param- 
eter, which requires £,s,norm < 1 for detection of entan- 
glement but in the presence of arbitrary number N . to 
show perfect overlay. Our proof justifies the normaliza- 
tion procedure used in recent experiments that report 
spin squeezing and particle entanglement 043 > f° r a re ~ 
pulsive BEC with fluctuating total numbers. 

The criteria Eq. (|28p and (|2l?|) enable an unambiguous 
detection of entanglement in the presence of number fluc- 
tuations. We note that a similar immunity of the Peres 
positive partial transpose PPT entanglement measure to 
loss was shown in Ref. 53], and for other entanglement 
measures to loss in Ref [54]. 

The figures include both attractive (g < 0) and re- 
pulsive (g > 0) interactions. In the attractive case, the 
criteria (|23H28p for entanglement are satisfied when ap- 
plied directly to the modes of the two wells, then de- 
scribed by a, b. In the repulsive case, an interferometric 
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Figure 2. (Color online) Entanglement of a ground state BEC 
including Poissonian atom number fluctuations. Here N = 100, for 
fixed gl k = 10 3 . (i) Entanglement between wells a and b is detected 
if Eph < 1 (blue and red dashed curves) or Ehz < 1 (equation 11231 
(purple and green solid curves). Curves a,b (purple solid and blue 
dashed) are for two- well BEC modes a,b; curves a,i n t,bi n t (green 
solid and red dashed) are for a, b formed from BEC modes a^, bi 
input to the M-Z interferometer sequence depicted in Figure 1. 
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If /it 

Figure 4. (Color online) Measured phase uncertainty A</> (1->2| i. 
for Poissonian fluctuations. Phase sensitivity is better than shot 
noise level if A<j> < 0.1 = l/y/~{N). The red dashed curve cor- 
responds to the critical value of (N)g/K ~ 43.6 which gives the 
lowest value of E v ^ for the repulsive regime. Increasing (N)g/n 
to 10 3 (dotted) and 10 4 (solid), improves the optimum phase un- 
certainty, but the useful region of phase angles become narrower. 
Optimum sensitivity is at tp ± tt/2 and is determined by the spin 
squeezing parameter §g norm . For detecting an unknown phase 
using only two orthogonal measurements, the sensitivity over the 
quiet quadrant ip/n = i ± j becomes relevant. 
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Figure 3. (Color online) Spin squeezing and phase measurement 
sensitivity parameters for ground state of a two-well BEC with 
(N) = 100. Plot of spin squeezing parameters fg, £5 (purple 
and green solid curves) and the normalized parameters £g h , ph 

(blue and red dashed curves), for Poissonian number fluctuations. 

Y 1 z y ' / ' z 

States with £ s ,£ s ' ^ < 1 show sub-shot noise enhanced phase 
sensitivity for measurements o f sma ll rotations about a fixed phase. 
Inset shows the asymptotic ((A</>)^/ 2 ~ 0(\/2/N)) behav- 

ior with increasing N for fixed g/n = 10 3 . The Heisenberg limit is 
~ 0(y/T/N) m. 



sequence is necessary: two- well modes aj,6j are phase- 
shifted and placed through a MZ beam-splitter. The 
reason for this is that while an attractive BEC reduces 
fluctuations directly in the plane of J X ' Y , a repulsive 
BEC reduces fluctuations in a different plane, namely in 
J X ' Z . Without the additional input phase-shifter, the 
phase-measurement is in the Y — Z plane, where only 
one phase has reduced fluctuations, as observed [5|- 



B. BEC interferometric phase measurement 

Now we turn to the question of using the BEC two- 
mode states for the purpose of interferometric phase mea- 
surement. The spin squeezing measure £s in the presence 
of fixed N, and the normalized parameter £,s,norm in the 
presence of fluctuating numbers, give sensitivity for mea- 
surements of small rotations about a fixed phase below 
the SQL when £ s < 1 and £ s ,„ orm < 1 (Eq. ([55])), The 
figure reveals this spin squeezing to be predicted for 
a wide range of parameters of the ground state solution. 
The inset shows the reduction in noise to be near Heisen- 
berg limited with the scaling ~ 0(^2/N) evident. 

The worse-case sensitivity of the interferometer to an 
arbitrary angle defined within the two quiet quadrants 
of measurement for an unknown phase cj> is given by the 
rj p h of Eq. (|34j) . which is minimized according to the 
phase entanglement measure E p h of figure 01 The best 
scaling of E p h with TV is given as J 2 / 3 , that of the Cj 
coefficients, and is achieved at the critical value of Ng/n. 
This implies, from the phase sensitivity measure r] p h, a 
sensitivity of (A<j)) worse ~ 0(N~ 2 / 3 ). 

For the proposed BEC interferometer, one can evaluate 
the actual range of sensitivities for the unknown incom- 
ing phase (f) using Eq (|32j) directly. The different range of 
phase-noise reduction as a function of measured phase- 
angle and BEC interaction strength is shown in figure [4] 
The best sensitivity is obtained at y> = tt/2 and the value 
for A(f> is determined by the spin squeezing parameter 
Eq. (|29p t;s,norm, which reduces to £5 for fixed number 
N. Where one measures an unknown phase cf> using only 
two orthogonal measurements (8 = and 8' = tt/2). only 
the sensitivity over the "quiet" quadrant indicated on the 
graph by the region (p/ir = 1/2 ± 1/4 becomes relevant. 
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In this case, the worse case sensitivities are at the edges 
(ip = tt/4 and (p = 37r/4) and determined by the value of 
r\ v h- : Eq. (|33| . This parameter is optimized by minimiz- 
ing the two-mode entanglement parameter E p h, which 
reduces to the Hillery Zubairy entanglement parameter 
when N is fixed. This demonstrates the trade-off between 
noise-reduction and range of measurable phase. 



VI. CONCLUSIONS 

In summary, we have introduced the normalized rel- 
ative phase quadrature operator J (</>) as the most di- 
rect operational expression of how interferometric mea- 
surements give rise to phase information. A correspond- 
ing phase-entanglement measure, as quantified by E p h, 
describes a useful physical resource for phase measure- 
ment. This directly quantifies the measurement sensitiv- 
ity increase above the standard quantum limit, as 
decreases towards a maximally entangled state. We also 
introduce a normalized phase squeezing measure, ^ p \, 
which signifies entanglement between qubits or particles. 
Both measures are normalized in terms of the total par- 
ticle number, using the generalized Moore-Penrose in- 
verse method. We show how it is possible to improve 
BEC phase measurements so that a range of unknown 
phases have reduced phase-noise, rather than just one 
pre-selected phase. 



APPENDIX 

In this Appendix, we provide a detailed proof of the 
normalized entanglement criterion, equation (|28[) and the 
normalized spin-squeezing criterion, 



Phase-entanglement criterion 

We wish to show that entanglement between modes a 
and b is confirmed via an entanglement phase criterion, 
equation (|2"5)) . First, we note the general result that if N 
commutes with an arbitrary hermitian operator O having 
eigenvalues j, we can introduce a limiting procedure to 
obtain a normalized mean value, (O) = (ON + ), where 
N + is the generalized Moore-Penrose [35| inverse of the 
number operator N, so that: 



(ON+) = limX; 



nj 



■P(n,j) 



^2n + jP(nJ). 



(36) 



Here P(n, j) is the probability for simultaneous outcomes 
n, j for N and O respectively, and the eigenvalues of the 
generalized inverse operator N + are n + = rtT 1 for n > 0, 



with n + = for n = 0. Hence, the expectation value for 
the ratio becomes: 

(6) = ^<o)„p„, 

n>0 

where (0) n = £\ P(j\n)j , P n = ZjPi^j), P(j\n) = 
P(n,j)/P n) and we define (0) n = (0} n n + . Similarly, 
the corresponding variances, A 2 = (O 2 ) — (O) 2 , can 
be expanded as: 



A 2 = J2(0 2 )nPn 



However, we know from elementary variance properties 

that Y. n (0) 2 n P n > [Z)„(0)«-Pnj , hence we can write 
that: 

n n 



o 



n + ) P n 



(37) 



Here, as usual, we have defined A 2 = J2j P(j\ n )j 2 ~ 
(0)1 

Next, we apply this result to normalized spin variances, 
giving: 



(A.F) 2 + (A,F) 2 > E P n [A 2 n J x + Air], 



where we have used the definitions that: 

A 2 ,F = E p ™( JX2 )™-( JX ) 2 - 



(38) 



Finally, if we assume the separability equation (|2"2")l . we 
see that for a separable density matrix with a fixed total 
particle number: 

A 2 J* + A 2 J^>n/2, 
which means that for normalized operators, 

(AF) 2 + (AJ>-) 2 >iEM- + ) 2p « 

n 

= ^E- +P » 

n 

Here we note that the generalized inverse has many 
of the properties of the standard inverse, in particular 

that N (^ + ) = N + , and that of course no singularities 
occur with this criterion. When this condition is violated, 
we must have an entangled state, which leads to 
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Phase-squeezing criterion 

Next, we wish to show that entanglement between N 
spin- 1/2 systems, where the number of systems can fluc- 
tuate, is confirmed by the normalized spin squeezing cri- 
terion of equation ([2"9"[) . For TV spin 1/2 separable sys- 
tems, where N is fixed and nonzero, one finds that [381 1 : 



(AJ z Y>^[(J x f 
> N+(J X ) 2 



> J^[(J X ) 2 + (J Y ) 2 } 



The last expression, using the generalized inverse of N, 
holds even when N = 0. Thus, using the normalized 
variance inequality (|37p . we obtain: 

A^>£P„(n+) 2 [A*n (39) 

n 

>Y,Pnn + (J x N + )l (40) 

n 

Therefore, using the result that n 2 n + = n, we se that: 



(JV)A 2 ?>(E^ V } 



n>0 



n>0 



Next, we use the Cauchy-Schwarz inequality: 
{E„>o a; n>o}{E„>o2/n} > IEn>o z n»n| 2 where 



%n = V n Pn and y n = y P n n + (j x N + ) n , and hence: 



(N)A A J Z > 



.71 >0 



(J X N+) =|(J 



rX\|2 



This proves the phase-squeezing criterion, equation ([29 
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